Microscopic Theory of a Quantum Hall Ising Nematic: Domain Walls and Disorder 
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We study the the interplay between spontaneously broken valley symmetry and spatial disorder in multivalley 
semiconductors in the quantum Hall regime. In cases where valleys have anisotropic electron dispersion a previ- 
ous long- wavelength analysis [Phys. Rev. B 82, 035428 (2010)] identified two new phases exhibiting the QHE. 
The first is the quantum Hall Ising nematic (QHIN), a phase with long-range orientational order manifested in 
macroscopic transport anisotropics. The second is the quantum Hall random-field paramagnet (QHRFPM) that 
emerges when the Ising ordering is disrupted by quenched disorder, characterized by a domain structure with a 
distinctive response to a valley symmetry-breaking strain field. Here we provide a more detailed microscopic 
analysis of the QHIN, which allows us to (i) estimate its Ising ordering temperature; (ii) study its domain- wall 
excitations, which play a central role in determining its properties; and (iii) analyze its response to quenched dis- 
order from impurity scattering, which gives an estimate for domain size in the descendant QHRFPM. Our results 
are directly applicable to AlAs heterostructures, although their qualitative aspects inform other ferromagnetic 
QH systems, such as Si(l 11) heterostructures and bilayer graphene with trigonal warping. 
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I. INTRODUCTION 

Two-dimensional electron gases (2DEGs) placed in high 
magnetic fields exhibit a multitude of phases associated with 
the quantum Hall effect (QHE).^ Particularly interesting in 
this context are situations in which quantum Hall ordering 
is accompanied by the breaking of internal symmetries - 
such as the global symmetries associated with the electron 
spin,Elor valle>^ "^ or layer^ pseudospin. The resulting broken- 
symmetry state, termed a quantum Hall ferromagnet, pos- 
sesses in addition to the topological order common to all quan- 
tum Hall states a distinctive set of phenomena relating to the 
low-energy pseudospin degrees of freedom. These include 
charged skyrmions, finite-temperature phase transitions, and 
Josephson-like effects, to name a few. 

Recent experimental^K^ and theoretical work has focused 
on the case in which the symmetry in question is between 
the different valleys (i.e., conduction band minima) of a semi- 
conductor. In previous worli^ involving two of the present 
authors, it was noted that a generic feature of such multival- 
ley systems is that the point-group symmetries act simulta- 
neously on the internal valley pseudospin index and on the 
spatial degrees of freedom. This linking of pseudospin and 
space has significant consequences at "ferromagnetic" filling 
factors, such as i/ = 1: 

(i) in the absence of disorder, pseudospin ferromagnetism 
onsets via an Ising-type finite-temperature transition and 
is necessarily accompanied by broken rotational symme- 
try, corresponding to nematic order. The resulting state at 
T = 0, dubbed the quantum Hall Ising nematic (QHIN), 
has an intrinsic resistive anisotropy for dissipative trans- 
port near the center of the corresponding quantum Hall 
plateau. 

(//) as a quenched random field is a relevant perturbation to 
Ising order in (i = 2, the QHIN is unstable to spatial dis- 
order - such as random potentials or strains - that gives 
rise to such fields. Disorder thus destroys the long-range 



nematic order, giving rise to a paramagnetic phase. Pro- 
vided that there is (arbitrarily weak) intervalley scatter- 
ing, this continues to exhibit the QHE at weak disorder 
and low temperatures, and is hence termed the quantum 
Hall random-field paramagnet (QHRFPM). Transport in 
this phase is dominated by excitations hosted by domain 
walls between different orientations of the nematic order 
parameter, and is extremely sensitive to the application 
of a symmetry-breaking 'valley Zeeman' field - for in- 
stance, due to uniaxial strain -which can tune between 
percolating and disconnected domain walls. 

Two aspects of this picture are particularly striking and 
should apply to a variety of valley quantum Hall ferromag- 
nets. The first is the role of valley anisotropy in establish- 
ing the nature of the symmetry breaking. Systems with val- 
leys that are isotropic (for instance, graphene), or have iden- 
tical anisotropics (such as Si (110) quantum wells) will ex- 
hibit an enhanced SU(2) valley pseudospin symmetry. It is 
the valley anisotropy in the present situation that entangles 
rotations in space with those in pseudospin space, and also 
reduces the order parameter to an Ising variable. Similar be- 
havior is expected for bilayer graphene once trigonal warp- 
ing of the band structure is included, and for Si (111) het- 
erostructures. Second, we emphasize that the QHIN and the 
QHRFPM that naturally emerge in this situation both exhibit 
quantum Hall behavior, but on parametrically different scales: 
the QHRFPM shows quantized conductivity only at temper- 
atures below the scale of domain wall-excitations, typically 
dominated by weak interactions and/or disorder, and hence, 
much lower than the intrinsic anisotropy scale characteristic 
of QH transport in the QHIN. 

A specific example of experimental interest^ ^ ^ ^^ and our 
focus in this paper is the case of wide quantum wells in AlAs 
heterostructures. Here, two valleys with ellipsoidal Fermi sur- 
faces are present, as shown in Fig. [T] Owing to the anisotropic 
effective mass tensor in the two valleys, individual electronic 
states no longer exhibit full rotational invariance. Only dis- 
crete rotations of the axes, accompanied by a simultaneous 
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FIG. 1. (a.) Model band structure used in this paper, appropriate 
to describing AlAs wide quantum wells, (b.) Different phases as 
determined by comparing Imry-Ma domain size ^/m to sample di- 
mensions Ls. Top: For ^/m ^ Ls we find the QHRFPM. Bottom: 
For ^iM ^ Ls the system is dominated by the properties of a single 
domain, and is better modeled as a QHIN. At intermediate scales, 
Ls "^ ^iM there is a crossover. 



amount of 'texturing' in a soft domain wall as a function of the 
anisotropy (Sec. IIIE| ) - we note that texturing is a prediction 
of the NLcrM that remains valid at high anisotropics. We next 
turn to an analysis of disorder within the microscopic theory, 
where we first establish that anisotropics in the screened ran- 
dom impurity potential act as a valley- selection mechanism, 
translating into a random field acting on the Ising order param- 



eter (Sec. IV A ), which we compute in Landau-level mixing 
perturbation theory. We discuss how to estimate the strength 
of the disorder from the mobility, a measure that is readily 



accessible to experiments (Sec. IVB ). Taken together, the do- 
main wall parameters and the random field studies yield es- 
timates for the characteristic domain size due to the disorder, 
allowing us to make partial contact with experiments (Sec. 
\V\ . All these results are obtained for the microscopies of the 
AlAs heterostructures which were the original motivation for 
our study of valley-nematic order. However, the qualitative 
features of domain wall structure, random-field disorder, and 
dipole moment physics apply mutatis mutandis to other mul- 
tivalley systems in the QH regime. 



interchange of the valleys remain as symmetries of the system. 
It is in this specific sense that the internal index is entangled 
with the spatial symmetries. 

The existence of the two phases was originally established 
within a long-wavelength nonlinear sigma model (NLaM) 
field theory, which also provides a caricature of their prop- 
erties and the above phase diagram in the weak-anisotropy 
limit. While it is expected that this treatment captures qualita- 
tive features of valley Ising physics reasonably well, to make 
a quantitative connection to experiments a microscopic un- 
derstanding is essential. Here, we provide such a microscopic 
analysis of the QHIN, focusing specifically on properties of 
domain walls which as we have argued are central to this sys- 
tem. 

A summary of the main results of this paper, which also 
serves to outline its organization, follows. We first place this 
work in context by providing a summary of the important as- 
pects of valley-nematic ordering in the quantum Hall effect 
in Sec. [ll| focusing on qualitative features of the phase di- 
agram, the role of thermal fluctuations and quenched disor- 
der and thermal fluctuations, and transport signatures of the 
QHIN/QHRFPM phases. We then proceed to our technical re- 



sults. First, we set up a Hartree-Fock formalism (Sec. Ill A), 
which we use to obtain a mean-field estimate of the transi- 
tion temperature out of the thermally disordered phase (Sec. 
IIIB| ). We proceed to construct a solution of the HF equations 



Iirc| ), where 



corresponding to a 'sharp' domain wall (Sec. 
the valley pseudospin changes its orientation abruptly at the 
wall; this is expected to be an accurate description of phys- 
ical domain boundaries in the 'strongly Ising' limit of large 
mass anisotropy. We determine the properties of the sharp 
wall as a function of the mass anisotropy, specifically its sur- 
face tension and dipole moment, the latter a property which 
is not captured in the NLaM limit. We clarify the effect of 
this dipole moment on critical behavior and domain wall en- 
ergetics (Sec. |IIID| ). We then relax the sharp- wall approxima- 
tion and numerically solve the HF equations to quantify the 



II. OVERVIEW: PHASES, TRANSITIONS, TRANSPORT 

The temperature-disorder phase diagram of multivalley 
2DEGs exhibiting Ising valley ordering can be sketched as fol- 
lows (see Fig. [2]). In the absence of disorder, there is a finite 
temperature transition into an Ising nematic ordered phase, 
which exhibits transport features of the QHE. While strictly 
speaking, the QHE is a zero-temperature phenomenon, in a 
slight abuse of terminology we will nevertheless refer to the 
entire phase below Tc in the zero-disorder limit as the QHIN. 
The quantization of the Hall conductivity and the vanishing of 
the longitudinal conductivity are only exponentially accurate 
at finite temperature. While there is a thermodynamic transi- 
tion associated with the Ising valley ordering, the conductiv- 
ity exhibits a crossover rather than a singularity at Tc. The 
orientational symmetry breaking of the Ising nematic phase 
is reflected in the anisotropic longitudinal conductivity of the 
QHIN where (Jxx/ayy ^ 1. Upon adding disorder, the Ising 
transition is destroyed and at T = the system is in the 
QHRFPM phase. Above this at finite temperature (shaded re- 
gion in Fig. [2]) we once again find exponentially vanishing lon- 
gitudinal conductivity and exponentially quantized Hall con- 
ductivity, but the response is now isotropic: dxx/o'yy = 1. 
With similar caveats as in the clean case we will refer to the 
entire shaded region above the T = line as the QHRFPM. In 
contrast to the QHIN, there is no thermodynamic phase tran- 
sition into the QHRFPM at T > 0, only a crossover in the 
conductivity at a temperature scale T* (dashed line in Fig. [2]) 

We emphasize that there is an important qualitative differ- 
ence between the QHIN and the QHRFPM, over and above 
the anisotropy in the former. Namely, the crossover into a 
quantized Hall response in transport is governed by differ- 
ent physical mechanisms. In the QHIN, this crossover occurs 
at a scale set by the exchange energy, effectively the single- 
particle gap. Asp ~ e^ /cIb in the QH ferromagnetic ground 
state. This also sets the scale of the Ising T^ upto a numeri- 
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FIG. 2. Phase diagram as function of temperature (T) and disorder 
strength (W), showing behavior of conductivity. The phases and 
critical points are defined in the introduction. 



cal factor that depends on the mass anisotropy. In contrast the 
QHRFPM is, as we have noted, characterized by multiple do- 
mains of differing Ising polarization. Here, the lowest-energy 
charged excitations are localized on one-dimensional domain 
boundaries, ^^ which in the strong-anisotropy limit can be un- 
derstood in terms of a pair of counterpropagating QH edge 
states of opposite pseudospin. The stability of the QHE then 
rests on the gap to creating domain- wall excitations. As this 
is induced by weak pseudospin symmetry-breaking terms in 
the Hamiltonian from both disorder and interactions, it is ex- 
pected to be small and the concomitant conductance quantiza- 
tion is thus fragile. At weak disorder, the dominant source ^^ 
of symmetry breaking is from intervalley Coulomb scattering, 
Viy , which thus sets the domain- wall gap Adw and hence the 
crossover scale T*. For sufficiently strong disorder above a 
critical strength W^ the energy gap stabilizing the QHRFPM 
collapses via the Fogler-Shklovskii scenaricP^ originally de- 
vised to describe the collapse of spin- splitting in quantum Hall 
ferromagnets in GaAs quantum wells. 

Whether a particular experimental sample will display 
the transport anisotropy characteristic of the QHIN, or the 
isotropic domain-wall dominated transport of the QHRFPM 
is a matter of quantitative detail, determined by the compara- 
tive energetics of the Ising exchange energy and the disorder. 
Their competition sets a characteristic 'Tmry-Ma' ^^ domain 
size ^iM in the random-field phase. The question then turns 
on whether the system consists of a single Ising domain or 
multiple domains, i.e. it depends on how the domain size 
compares to the sample dimensions, Ls (see Fig. [T]). The 
exchange strength is determined by the electron-electron in- 
teractions, while for the heterostructures of interest the disor- 
der is sensitive to the density of dopant impurities and their 
typical distance from the plane of the 2DEG. The effective 
mass anisotropy is important to estimates of both these quanti- 
ties, for in the isotropic limit there is a full SU(2) pseudospin 
symmetry, and potential disorder does not exhibit a prefer- 
ence for any particular pseudospin orientation. Thus, accurate 
estimates of these quantities picture are essential to make a 



FIG. 3. Valley symmetry-breaking field permits transport to probe 
the energy scales of the QHIN/QHRFPM. (Inset) Domain structure 
as function of disorder strength and valley splitting; dashed line 
shows a representative path in A^ leading to a transport signature 
similar to that in the main figure. AJ is the valley splitting for which 
the system is single-domain dominated. 



quantitative connection with experiments. 

The introduction of an externally applied valley Zeeman 
field - experimentally achieved via application of uniaxial 
strain to the 2DEG - provides a convenient probe of the trans- 
port scales in the QHIN and QHRFPM. First, this field intro- 
duces a single-particle splitting between valleys A^ that sta- 
bilizes the Ising nematic against the effects of disorder. Thus 
for sufficiently weak disorder and sufficently large A^, the 
anisotropic longitudinal conductivity should be clearly estab- 
lished. Second, in the case when for A^ = the disorder 
is sufficient that the sample is in the QHRFPM with multi- 
ple domains (for instance, along the dotted line in the inset 
of Fig. [Sj, application of the valley Zeeman field causes a 
crossover in the the longitudinal conductivity as a function of 
A^. A sketch of this is provided in Fig.|3j and can be under- 
stood as follows. For a disorder strength corresponding to the 
dotted line in the inset, the system crosses over from multi- 
ple domain to single domain behavior. This is reflected in the 
activation gap for longitudinal transport: in the multiple do- 
main regime, the gap is dominated by the domain wall scale 
/S^dw Deep in the single-domain regime, the gap is essentially 
set by the single-particle gap, which itself scales linearly with 
A^ ; the intercept of the asymptotic linear dependence can be 
used to extract the characteristic single-particle energy scale 
at zero Zeeman splitting. The sharp crossover between the 
two regimes can be understood qualitatively in terms of tun- 
ing domain walls in a random-field Ising model away from 
percolation by applying a constant symmetry -breaking field. 



III. MICROSCOPIC THEORY 

We will begin by developing a microscopic theory of the 
QHIN using the Hartree-Fock (HF) approximation.^^ We will 



focus on the case relevant to AlAs, with two valleys denoted 
by index /^ = 1, 2 and centered at Ki = {Kq^ 0) and K2 = 
(0, ifo) respectively, with mass anisotropyA^ = ^^ = ^^^ 
(a schematic dispersion is sketched in Fig.[T]) In each valley, 
the single-particle kinetic energy is 
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Working in Landau gauge, A 
lowest LL eigenfunctions are 



(0, -Bx), we find that the 
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1/4 



u^(:i^-Xy 



(2) 



Here, ui = l/u2 = A, and we have labeled states within 
a LL by their momentum py, which translates into a guiding- 
center coordinate via X = Pyi%. Henceforth, we shall ac- 
count for the spatial structure of the LL eigenstates by the 
standard procedure of projecting the density operators onto 
the lowest LLP 



A. Hartree-Fock Formalism 

We consider a rectangular system of dimensions L^^Ly. 
Since we are interested in sl u = 1 state, the total number of 
electrons in the system (which we take to be even for conve- 
nience) is N = N^ — ^^^2^ . For periodic boundary condi- 



27r£2 



tions in the ^/-direction, the guiding center coordinate along 
the x-direction is given by X^ 



27r^ 



between 



' 2 



1 and 



N 



^ n, with n an integer 
For the sake of brevity, we shall 

I 



continue to label states by X, but with the understanding that 
it is now a discrete index. Unless otherwise mentioned, all 
sums and products are over the full (finite) range of X. 

In addition to the electron-electron interaction energy, the 
lowest Landau level Hamiltonian must include the energy of 
the electrons interacting with the potential of the positively 
charged background, which depends on the form of the back- 
ground charge density. We shall take the positive charges to 
have orbitals of the same form as electronic states in the two 



valleys, and corresponding occupation numbers nk 



p,{v) = Y,n^'^rx,A^)^xA^) 



ib). 



X,K 



(3) 



1 



.W 



1. Note that for 



with n = 1, 2 as before and n 

any choice of ni^ ^ satisfying the latter constraint, ^^(r) is the 
same uniform constant. However, a judicious choice of the 
background charges will allow us to cancel divergences of the 
Hartree contribution, as we will see below. 

In order to model boundaries between valley domains we 
also add a spatially varying single-particle pseudospin split- 
ting that increases linearly in X from negative to positive 
across the system^^ which models the external random valley 
Zeeman field from disorder. This serves a twofold purpose: 
first, it pins the domain wall^^ near n = 0, which is desirable 
for a stable numerical solution even in the clean limit; second, 
it allows us to study how the domain wall properties change 
as we vary the characteristic length scale and typical strength 
of the random field that leads to domain formation. 

With these preliminaries, the second-quantized Hamilto- 
nian projected to the lowest Landau level can now be written 
in the Landau basis: 



^ - 2 Z^ Z^ ^i^'X',kX^kY^k'Y'^i^'X'(^i^X Z^Z^Y^l ^IX.KY^KY^'^y ^^2 ^2X,kY ^kY^i^Y 
k,k' X,Y k X,Y 

X' ,Y' 



(4) 



The first term is the electron-electron interaction, the second 
is the interaction between the electrons and the positive back- 
ground, and the third term is the single-particle splitting. As 
discussed in the next section, the characteristic energy scale 
of this is A^^ , and it varies over a characteristic distance d 
corresponding to the correlation length of the random field; 
rewriting this carefully, leads to the expression given, with 

- (which has units of energy). Note that because 
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2d 



d ^ Ib^ this term is fairly small even at the two ends of the 
system, where it is maximal. The final term is the self-energy 
of the background charge distribution, a positive constant that 
we omit forthwith. In writing Q, we have ignored 'umk- 



lapp' terms that lead to a net transfer of electrons between 
valleys (as these are exponentially suppressed in a/ls^ as well 
as terms that exchange a pair of electrons in the two valleys 
(suppressed by a factor of (a/^B)^)- At the scales of interest, 
even the latter term only contributes a small energy correction 
(^ 1% of the terms kept), and we do not expect their inclusion 
to significantly alter our conclusions. 

The matrix elements of the Coulomb interaction are given 
by the usual second-quantized form: 



ttKX^k'X' 
^K,'Y',K,Y 



d^rd^r'r.xir)r.'X'{r')V{r- 



xtp^,Y'{r')ip^Y{r) 



(5) 



where the single-particle wave functions were defined in the 
previous section. Note that momentum conservation requires 

that X + X' = r + r^ 



1. Energy scales 

Throughout the remainder of this paper, we present our re- 
sults in dimensionless units. We measure energy in units of 
the Coulomb energy e^ /cIb where e is the dielectric constant 
appropriate to the hetero structure under consideration. For the 
AlAs devices which are our primary focus, e ~ 10. The mag- 
netic length, Ib ^ 8nm for a magnetic field of lOT, so that 
e^ /els ^ 200 K. The surface tension of Ising domain walls is 
measured in units of e^/e£^, roughly 25 K/nm for this choice 
of parameters. 




FIG. 4. Mean-field and NLcrM estimates of Tc. Dashed line shows 
the anisotropy (A^ ?^ 5.5) appropriate to AlAs. 



B. Estimates of Tc 

Our first application of the microscopic theory will be to es- 
timate the Ising ordering temperature Tc for the clean system 
via finite-temperature Hartree-Fock theory. A standard mean- 
field decoupling of the Hamiltonian ^ in the density channel, 
(c^j^c^/y) = n,^6,^,^'5xY^ where the occupation numbers are 
assumed independent of position, yields 
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MF 



~9 2^ Z^ ^^^^Y^^X <Y^^y 



X t 
Y^i^Y 
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(6) 
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We simplify the Hartree term by taking n)^^ = n^. For N -^ 
oo, we may use translation invariance of the potential to write 

^MF = Zlx K ^K.clx(^K.x, where 



1 Y^ TrKX,KY 

k'Y 



2^' 

Y 






is independent of X. 

We seek a solution where the ground state spontaneously 
breaks valley symmetry; without loss of generality we may 
assume it is polarized in valley 1, and take the energy splitting 
to be A, whence 



ni 



gA/fesT 
l^gA/Zc^T' 



1 



^2 



oA//cbT 



(8) 



For this ansatz the self-consistency condition corresponds to 
A = e2 — ei = 74ini — A2n2, where after a tedious calculation 
we find (defining 1 = 2,2 = 1) 



A, 



ttkX^kY 
~^K,Y,RX 



-T 

9 2^ ^kX,kY 
^ Y 



TrK,X,K,Y 

■ ^kX,kY 



V 



kX,kY 
kY,kX 



tkX.kY 



(9) 



(This cancellation of the Hartree contributions from the two 
valleys is the reason for the choice of background charge made 



previously.) Valley symmetry requires that Ai = A2, which 
yields the self-consistency condition A = A^ tanh 2a~t • ^^ 
the standard comparison of the slope of both sides of this 
equation at A = 0, we find for the (mean-field) transition 
temperature 



kBT^' = \a. 



4 Z^ ^kX.kY 
Y 



1 e" 



dx 



dy- 



-h(^^^y') 



1 e^ K{1- l/A^^ 



2(2^)3/2 eiB VX 



(10) 



where K is the complete elliptic integral of the first kind.!^ 
Note that this mean-field expression for Tc has some un- 
physical aspects - most notably it is nonzero even in the 
Heisenberg limit (A -^ 1), and decreases with increasing 
mass anisotropy. This will be corrected in an RPA spin- 
wave calculation of quadratic fluctuations about the mean- 
field ground state. In particular, the fluctuations drive T^^ to 
zero in the isotropic Heisenberg limit. Furthermore, as spin- 
wave gap scales roughly with the Ising anisotropy, the debil- 
itating effect of spin waves on T^^ is suppressed at strong 
anisotropy, offsetting the decrease in the energy scale pre- 
dicted by the mean-field theory. As the spin-wave calcula- 
tion is technically involved and not too informative, we pro- 
vide instead an alternative estimate of Tc for comparison: 
T^ ~ Airps log"^[ps/a^^], obtained from the NLaM with 



stiffness ps 



0.025- 



and Ising anisotropy a 



16V27reiB 

0.01 -73- (^ ~ 1)^' whose leading dependence of A was com- 
puted in a gradient expansion in Ref . [12] We plot both esti- 
mates inFig.|4] 



C. Properties of Sharp Domain Walls 

We turn now to an analysis of 'sharp' domain walls. These 
are solutions to the HF equations where the valley pseudospin 



abruptly changes orientation from one Landau gauge orbital 
to the next. We will determine the properties of the sharp do- 
main wall as a function of the anisotropy. While analytically 
tractable, this approximation is expected to be a good descrip- 
tion of the domain wall only at strong anisotropy, but never- 
theless provides a valuable complementary perspective of its 
properties in a regime where the NLaM is no longer valid. If 
we take as the ground state a fully pseudospin polarized Slater 
determinant with all the electrons in valley 1 : 



\^g)=1[^\x\0) 



(11) 



X 



then a domain wall is captured by a Slater determinant of the 
form 

I^Dw) = n (^^4x + ^xclx) |0). (12) 

X 

The sharp wall corresponds to the case ux = IjVx = for 
X < and ux = O^Vx = 1 for X >. We once again consider 
the Hamiltonian ^ with ^ = and assume a background 



charge distribution polarized in valley 1, i.e. 



.(^) - 



1. Two 



properties of the domain wall will be of especial interest to us: 
its dipole moment and its surface tension. 




FIG. 5. (a.) Surface tension and (b.) dipole moment of a sharp DW 
as a function of the effective mass anisotropy. Dashed line shows the 
anisotropy (A^ ?^ 5.5) appropriate to AlAs. 



1. Surface Tension 

The first quantity of interest is the domain wall surface ten- 
sion - the energy per unit length of the wall. This provides 
a measure of the Ising exchange energy appropriate to the 
strong-anisotropy limit. Note that within the NLaM the do- 
main wall surface tension depends both on the stiffness and 
the Ising anisotropy. In the microscopic theory, we find the 
surface tension (energy per unit length along the wall) of a 
sharp domain wall to be the sum of three contributions: 



a(A)= lim ^((i/)DW 



Eo) 



= El -\- Ell 
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(13) 



Here, Eq and {{H)j)w — ^o) are the energies of the ground 
state and the sharp domain wall. The three contributions are 
individually convergent, and can be written as follows. The 
first term. 
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measures the Hartree cost, and can be simplified as 



El 



where 



1 e^ r . . . (^ . fx{x)fx{x') 
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(15) 



fx{x) = erfc (-xA-^/2) - erfc (-xX^^^) (16) 



with erfc the complementary error function. The second term, 
^n = - Yl ('^ix,ix' - ^2X,2X' J (1'7) 



X,X' = 1 



is the difference in the 'bulk' exchange energy between 
ground state and domain wall state from orbitals near the cen- 
ter or the edge.The final contribution measures the loss of ex- 
change energy since the two valleys have vanishing exchange 
matrix elements: 



E^ 



III 



E E^. 



IX, IX' 
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(18) 



We find cr( A) by numerically computing the convergent sums 
£^i, Ell and ^ni, in each of which we can truly take the upper 
bounds on X to infinity. Note that a (A) depends logarith- 
mically on Ly from the upper bound in the integral in (15). 
Ignoring this weak dependence, we can take Ly ^ oo in 
integrals over Qy and numerically integrate each term to ob- 
tain the surface tension as a function of anisotropy, plotted in 
Fig.[5](a). 



2. Dipole Moment 

Consider for a moment a long-wavelength description of 
an Ising nematic in terms of a single-component Ising order 
parameter field (p, and consider a domain wall parallel to the 
y-8ixis at X = 0, between regions with opposite Ising polar- 
ization (i.e., (p ^ ±1 8iS X ^ ^oo.) The remaining rotational 
symmetry is a rotation R^^ that takes x -^ —x^y -^ —y. For 



the given configuration, we have Rt^(p{x, y) = (/^(— :r, —y) = 
—(^{x^ y). Observe that under this symmetry, dx(p{x, y) is left 
invariant. From this it is not too difficult to show that V(/?(r) 
transforms as a vector under a rotation by tt, and thus has 
the same symmetry as that of a dipole moment normal to the 
domain boundary and oriented in the direction of decreasing 
Ising polarization. It is quite straightforward to find a micro- 
scopic origin for the dipole moment. Recall that the spatial 
extent of the Landau gauge orbitals in the X-direction is dif- 
ferent in the two valleys. At a domain wall, the charge distri- 
bution from valley 1 decays with a smaller Gaussian envelope 
than the growth of charge from valley 2. Assuming a uniform 
positive background, this leads to a dipole moment associ- 
ated with the interface between the two valleys and oriented 
as above. Therefore the theory of an Ising nematic should 
properly include long-range interactions between dipoles tied 
to gradients in the Ising order parameter. However, these ap- 
pear only at higher orders in the gradient expansion than those 
used to obtain the leading terms in the long-wavelength the- 
ory and represent a small perturbation in the weak-anisotropy 
limit. It is easier to compute the dipole moment at a domain 
wall within the microscopic theory: for our choice of back- 
ground charge, it is straightforward to show that the charge 
distribution associated with a sharp domain wall is 

Aot(r) =Pe(r)+Pbg(r) 

N/2 

= - E {("^ - l)^ix(r)^ix(r) 

X=-Ar/2+l 
N/2 

= - E Mx(r)^ix(r) - ^2*;,(r)^2x(r)} . 



X=l 



(19) 



Performing the summations and using the explicit form of the 
single-particle wavef unctions, we can verify that Ptot(r) cor- 
responds to a pair of dipolar charge distributions, one located 
at the domain wall {X = 0) and the other at the right edge 
of the system (since the background falls off with a different 
exponential than the electronic density.) Some care must be 
taken to separate the contribution of just the dipole at the cen- 
ter so that we have a controlled L^^ ^ oo limit; after some 
work, we find the dipole moment per unit length of a domain 
wall is given by 






1 r 



t^e~^ dt 



(20) 



Note that the dipole moment changes sign under A -^ 1/A, 
reflecting the fact its existence is directly tied to the mass 
anisotropy; we plot this in Fig. [5] (b). We reiterate that the 
dipole moment associated with the DW is a generic feature 
of an Ising model in which the two phases are distinguished 
by an orientational symmetry-breaking order parameter; how- 
ever it is not captured by the NLaM description of Ref. 12 at 
leading order in the limit of weak anisotropy. 



D. Does the Dipole Moment Matter? 

A central result of our microscopic study is that there is 
indeed a nonzero dipole moment at the domain wall as sug- 
gested by the symmetries of the system. However, as we have 
emphasized this physics is invisible in the weak-anisotropy 
NLcrM treatment on the basis of which we sketched the phase 
diagram of the system with temperature and disorder and dis- 
cussed qualitative features of these phases. As a consequence 
of this dipole moment, there are long-range interactions be- 
tween different portions of a domain wall and between dif- 
ferent domain walls. Do these perturbations to the original 
long- wavelength theory affect the physics? We will address 
two separate questions: the role they play at the Ising tran- 
sition in the absence of disorder, as well as the interplay of 
the long-range couplings with the formation of domains in 
the Ising phase. As both questions should have universal an- 
swers independent of the microscopic model, it will suffice to 
consider the role of the dipole-dipole interactions in the long- 
wavelength theory. Therefore we consider the free energy of 
the 2D Ising model. 
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and add to it a perturbation appropriate to a long-range inter- 
action between dipoles: 
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and determine its effect on the critical theory and domain for- 
mation with disorder. 

(i) Irrelevance at Tc. Using the fact that p^. ~ \/(p{r), we 
have for dimensional purposes 
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where we have ignored angular factors as we are really 
only interested in power-counting. Recall^^ that a long- 
ranged spin- spin interaction scaling as l/x^~^^ is irrele- 
vant at the short-ranged Ising critical point if a > 2 — t^sr 
where t^sr is the anomalous dimension of the Ising field 
in the short-ranged theory. For dipolar interactions in the 
d = 2 Ising model we have t^sr = 1/4 and a = 3, 
and thus ([22]) represents an irrelevant perturbation at the 
finite-temperature Ising critical point Tc. 

(ii) Imry-Ma domain formation at T = 0. Recall that the 
standard Harris criterion^^/Imry-McP'^ll argument in the 
2D Ising ordered phase proceeds as follows: we flip 
spins to orient with the random field to gain an energy 
(X L, at the cost of a introducing a smooth domain wall 
whose energy also scales as L; thus, for a sufficiently 
weak random field there is no advantage to introduc- 
ing domains. However, a more sophisticated argumenP^l 
notes that domain wall roughening can increase the en- 
ergy gain from the random field so that it scales as 



LlogL. Thus, disorder always destroys the Ising or- 
dered phase in d = 2. We have verified that long-range 
dipolar interactions do not affect the qualitative features 
of this argument, so that disorder remains a relevant per- 
turbation that destroys Ising order at zero temperature. 

Although the universal physics and the critical points are un- 
affected, one physical manifestation of the dipolar interactions 
is to increase the numerical value of the surface tension and 
thus renormalize the Ising stiffness upwards. As a conse- 
quence, the characteristic size of an Ising domain in the ne- 
matic phase is enhanced - note that owing to the exponential 
dependence of the domain size on the stiffness this can be a 
quite significant effect. 



E. Domain Wall Texturing 

Thus far we have focused on a sharp domain wall. Within 
the NLcrM, we find that domain walls are always textured: 
there is a length scale, set by the competition between the 
Ising anisotropy (that breaks the SU{2) symmetry down to 



Z2) and the stiffness. Does the texturing persist even when 
the NLcrM is no longer valid? We answer this partially via 
a self-consistent numerical solution of a domain wall, which 
reveals that some texturing does indeed persist into the strong 
anisotropy regime; we also study the texturing as a function of 
the random field gradient at the wall, as it provides additional 
information about how the domain wall structure is altered in 
the presence of disorder. 

We take ( pT) as the ground state as before, and the domain 
wall solution is given by ([12]) subject now to the constraint 
l^x P + l^x P = 1, and with the boundary condition that ux 
and vx approach 1 for X = {-N/2 + l)27ri% /Ly and X = 
N/2 X 2t:1\ /Ly respectively. This corresponds to a domain 
wall where the pseudospin rotates from valley 1 to valley 2 as 
we move from left to right. Note that, unlike in the sharp case, 
the wall is allowed to 'texture', i.e. cross over from one valley 
to the other over a finite length scale. In our simulations, we 
will take L^ = lOiriB^Ly = 30£b corresponding to A^ = 
150, and once again take the background to be fully polarized 
in valley 1 . 

Using Wick's theorem and the HF trial wavefunction in Q, 
we find 



f/«(X) + C/f(X)+5(^-i 




where the Hartree-Fock potentials are 

U^X) = E [VSA"^ (kxf - 1) + T/S'x^l.x'l^] , Uf{X) = Y: [^.S^ {\ux'\' - 1) + l^S'x^bx' 

X' Y 



(24) 
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In the above expressions we have subtracted off the energy of 
the ground state, so that we may consistently compare domain 
wall energies for different values of the anisotropy. 



IV. DISORDER IN THE MICROSCOPIC THEORY 



The optimization procedure proceeds iteratively, as follows. 
We begin with a trial wavefunction satisfying the boundary 
conditions, and in each iteration find the values of Um^Vm 
which optimize the HF energy, which are then used to gen- 
erate the HF potentials for the next iteration. Eventually, the 
procedure converges to a self-consistent solution. 



We estimate the degree of texturing by computing the mag- 
nitude of the x-component of the pseudospin in the domain 
wall configuration, since this is nonzero near the wall and van- 
ishes far from it. In Fig. |6J we plot contours of constant {Sx) 
in the anisotropy-field gradient plane, as well as the degree 
of texturing as a function of field gradient at A^ ^ 5.5, the 
anisotropy appropriate to AlAs. 



As discussed previously, disorder plays a central role in 
destabilizing the QHIN towards the QHRFPM. There are two 
primary sources of disorder: (i) random strains in the sys- 
tem can lead to a position-dependent shift of the energies in 
the two valleys- while the average strain (pseudomagnetic 
field) can be externally controlled, fluctuations of the strain 
are inevitable; and (ii) random fluctuations of the smooth elec- 
tric potential, Ud that arises from the screening of the po- 
tential due randomly placed donor impurities by electrons in 
the 2DEG also give rise to a random valley field. The ran- 
dom valley Zeeman field from the strain is difficult to quan- 
tify precisely, but is related to the anisotropy of the displace- 
ment field u{r) of the crystal from its equilibrium position: 
A^^^(r) oc {dx — dy)u{Y). The random electric field mech- 
anism can be understood via a straightforward application of 
perturbation theory and its value estimated from the sample 
mobility, as we now describe. 




FIG. 6. Domain- wall texturing from Hartee-Fock Theory. (Top) 
Contour plot of the average in-plane valley pseudospin (Sx) per unit 
magnetic length along the domain wall, as a function of the mass 
anisotropy A^ and the valley Zeeman field gradient g, with the latter 
on a logarithmic scale. The dashed line marks the anisotropy A^ ^ 
5.5 relevant to AlAs; note that there is still some texturing in this 
limit. (Bottom) Cut along dashed line, with ^ on a linear scale. 



A. Random Fields from Impurity Potential Scattering 

We briefly summarize the argument that leads to a coupling 
between a local anisotropy in the disorder potential and the 
Ising order parameter. Since the form factors of the two val- 
leys are different, we expect that the portion of the disorder 
potential that is antisymmetric in valley indices will lead to a 
spatially dependent single-particle splitting between valleys; 
in the limit when the cyclotron gap diverges, i.e. when the 
lowest Landau level approximation is exact, this is the only 
contribution, and we can argue from symmetry that the corre- 
sponding random field should take the form {d^ — dy)Ud (at 
least in the small-anisotropy limit). Note, however, that this 
term is a total derivative, and contributes significantly only at 
the boundary of a domain. To go beyond this, we must relax 
the cjc ^ oo limit, and allow for the effects of Landau-level 
mixing to first order in Ud\ since this allows for terms of or- 
der U^/huJc, the random field now receives contributions of 
the form {{dxUd)'^ — {dyUd)'^)/huJc, which is not simply a 
boundary term. 

To derive the higher-order contribution to the single- 
particle valley splitting from the Landau level mixing terms, 
we make a simplifying assumption: namely, we ignore in- 
teractions while computing the effect of mixing. While the 
interactions may combine with the effects of disorder to mod- 
ify details of the calculation, we expect that their neglect does 
not change the qualitative features of our results. The Hamil- 



tonian for noninteracting electrons in AlAs is, in the Landau 
basis 



^n 



n,X 



^ t/r(-q)e^^^''4,.,x^^n,.,x_,(25) 
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X ± 
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Here, we 



where we have defined X± — ^^ ^ ^ 
have expanded the notation of Section I to include Landau 
level indices n, m. In this basis, we have defined the ma- 
trix elements of the disorder potential Ud via U^'^{—c\) = 
/7rf(— q)F™(q), which naturally introduces the form factors 



rpnm 



(q) 




(26) 



for n > m, with F^^ 



(q) = F--(-q)*, where L" is the 
generalized Laguerre polynomial. 

Next, we compute a renormalized effective potentiaP^ 
within the lowest Landau level (where m = n = 0) by in- 
cluding Landau level mixing in perturbation theory. We find 

Ulll = E t^^eff (-q)e"^^''4.«.x,co...x_ (27) 

(l,X ,K 

where, to first order in Landau level mixing, 
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(28) 



We are primarily interested in the valley symmetry- 
breaking contribution from this term, so we consider only the 
portion antisymmetric in K:. Assuming that the disorder poten- 
tial is smooth on the scale of £b, we may expand in gradients 
of Ud', to quadratic order in QxjQy, only the n = 1 term in the 
sum contributes, and we find 

/7f (-q) = U'dM-q) - /7|,ff(-q) 

= -\{>^->^-')m-d'y)Ud]_^ 

+ ^ (A - X-') [{dxUdf - {dyUdf] ^ 

(29) 

The leading piece vanishes except on domain boundaries, as 
discussed; thus, the dominant valley splitting arising from im- 
purities is due to the second term. 

We focus our attention on a domain boundary, and as- 
sume that the distance between the centers of two domains is 
roughly the correlation length d of Ud- In this case, we simply 
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assume that the single-particle energy splitting changes sign 
linearly over a distance d, corresponding to the final term in 
Cb, with the overall energy scale A^^ set by the characteristic 
scale of the spatially varying random potential /7|^(r). 



B. Estimating Disorder Strength from Sample Mobility 

We may estimate the strength Ud of the smooth random 
potential from the measured sample mobility /i and the dis- 
tance d of the dopant atoms from the plane of the 2DEG, and 
using the results of the previous section, deduce the parame- 
ters of the random Zeeman field h. Taking the dopants to be 
Pois son-distributed, and assuming that the potential fluctua- 
tions are screened by electrons in the 2DEG, we can estimate 
the fluctuations of the potential^^ in the plane of the 2DEG to 
be 



\UM)\') = {Uodfe-''' 



(30) 



where Uq is determined by the screening length and should be 
proportional to the impurity density. 

The scattering rate due to this potential is 

977- 

A straightforward Boltzmann transport calculation of the 
transport relaxation time, assuming that it is dominated by the 
Fermi surface yields 



1 



(32) 
where the (1 — cosO) factor suppresses the contribution of 
small-angle scattering, which does not contribute to charge 
relaxation. For kpd^ 1, we have 



1 m 

Ttr TT?l 



^{Uodf- 



V^h^ 



8{mvFd)^' 
Using the fact that 1/rtr = e/{mfi) where /i is the mobility, 



(33) 



Uo 






(34) 



where we take m = ^Jrn^m^. 

Finally, we note that the characteristic length scale of the 
disorder potential is roughly the distance of the dopant plane 
from the 2DEG, allowing us to estimate that \sJUd\ ^ U^/d. 
Using the results of the previous section, the characteristic 
value of the symmetry breaking term A^^ is given by 



Af 



27rfi- 



r(m^ 



my)'^Ul 



(35) 



This result was used in Ref. 12 where it corresponds to a ran- 
dom field h ~ A|^/^^ in the NLaM. Since h is correlated 
roughly over a distance d, we find that the characteristic width 
of the random field distribution is I^ ~ {hdY\ this is the pa- 
rameter that quantifies the strength of disorder in our model. 



V. EXPERIMENTS 

As promised, we now turn to a discussion of probes of 
valley-nematic ordering via transport measurements. We will 
assume the ability to apply a valley- symmetry -breaking strain. 
Furthermore, we shall also assume that the maximal valley 
splitting that can be thus produced is sufficient to fully polar- 
ize the system in one of the valleys. We note that this is al- 
ready feasible for the samples studied experimentally thus far. 
We will also assume that the sample is engineered in a Hall bar 
geometry with principal axes parallel to the sample bound- 
aries, so that we may assume that the nematic anisotropy is 
oriented along the x- or y- direction of the sample. This re- 
moves ambiguity in the definition of components of the con- 
ductivity, but more importantly ensures that the anisotropics 
are observable in the Hall bar geometry. ^^ 

The cleanest probe of the valley ordering is to examine the 
longitudinal conductivity for anisotropy. A proxy for the ori- 
entational symmetry-breaking order parameter is the quantity 
C^ = (Jxx/cryy — 1. Note that it is important that both cFxx^ ^yy 
are measured in simultaneously, which can be conveniently 
accomplished in a four- terminal geometry.. The behavior of C, 
will exhibit quite distinct behavior as a function of tempera- 
ture and disorder strength, and will be affected by the applica- 
tion of a strain field. The principal distinction due to disorder 
is between 'clean' samples dominated by the properties of a 
single Imry-Ma domain, and 'dirty' ones which contain sev- 
eral domains. We identify four different cases: 

(i) Clean Sample, Zero Strain. Here, we expect that at high 
temperatures, the system is in the Ising thermal paramag- 
net phase, with no anisotropy, so C = 0; furthermore, C 
remains flat as the filling is tuned across the Hall plateau. 
As the temperature is lowered below the Ising Tc, the 
sample should enter the valley-ordered phase. Here, (" 
remains pinned to zero exactly at i^ = 1 i.e., the center of 
the Hall plateau. However, upon tuning the filling about 
u = 1, ( will change sign. This follows from the fact 
that the longitudinal conductivity goes from being dom- 
inated by hopping between hole-like levels of one valley 
to that between electron-like states of the opposite valley 
as the doping level crosses the center of the Hall plateau. 
The resulting longitudinal conductivities inherit the local 
anisotropy of Landau orbitals of the two valleys. ^^ ^^ The 
maximum value attained for T <^Tc can be estimated as 
Cmax ~ \\/mx/my - 1| = A- 1. 

(ii) Clean Sample, Under Strain. Application of strain to a 
clean sample should have little effect on the transport be- 
low Tc for one orientation of the strain, but should sup- 
press the anisotropy for the opposite orientation. In the 
paramagnetic phase, a strong valley polarization should 
result in transport signatures similar to that of the Ising 
ordered phase. 

(Hi) Dirty Sample, Zero Strain. For dirty samples, the 
anisotropy from the different domains cancel and we 
have C = for zero strain, at all temperatures. 
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(iv) Dirty Sample, Under Strain. Once again, application 
of strain to a dirty sample should polarize the system, 
and lead to transport signatures similar to the clean 
limit at zero strain, below Tc. As discussed previ- 
ously, the activation gap measured via longitudinal trans- 
port will be highly sensitive to the application of strain, 
and increase dramatically as the sample crosses over 
from multiple-domain to single-domain behavior and 
thus from domain-wall dominated to single-particle lon- 
gitudinal transport. 

As noted in the preceding section, the Imry-Ma domain size is 
exponentially sensitive to changes in microscopic parameters 
and thus estimating the domain size is a challenge. This can 
be circumvented to some degree by studying transport in sam- 
ples of different sizes and/or doping levels. For a given doping 
level, smaller samples are more likely to be in the clean limit 
as defined above, while lowering the doping level for sam- 
ples of a fixed size should weaken disorder to some extent. 
Also, the identification of clean and dirty samples is some- 
what loose; samples of intermediate size may show significant 
anisotropy even though there is no net Ising ordering, since the 
anisotropics of different domains may not fully cancel. 

Note that while the four-terminal probes are particularly un- 
ambiguous and striking, there is also useful information that 
can be gleaned from two-terminal transport measurements 
which only have access to a single longitudinal transport co- 
efficient. Here, the nematic symmetry breaking is encoded in 
the behavior of pxx as a function of the doping level. This 
will be minimal in the center of the Hall plateau, and grow 
as the filling is detuned from i/ = 1 in either direction. The 
mismatch in pxx for z^ < 1 and v > 1 will exhibit behavior 
similar to that described for C, in the different cases above. 

Finally, we note that random field Ising order is typically 
accompanied by a host of hysteretic effects^^ that might also 
be observable in experiments, particularly with an applied val- 
ley Zeeman field. 



VI. CONCLUDING REMARKS 

We have spent the majority of this paper focussing on a spe- 
cific instance of valley ordering relevant to experiments: the 
Ising-nematic order in AlAs quantum wells. As the reader no 
doubt appreciates by now, much of the richness of the phe- 
nomena discussed above stems from the inequivalence of the 
low-energy electronic dispersion in the two valleys. More 
specifically, the key observation underpinning our analysis is 
that the inequivalence between valleys is encoded by the fact 
that rotating between them necessarily requires a simultane- 
ous interchange of spatial axes; this has three striking conse- 
quences. First, in the presence of interactions the naive SU{2) 
symmetry associated with a generic 'internal' index is reduced 
to an Ising symmetry. Second, the intertwining of pseudospin 
and spatial rotations results in the transmutation of quenched 
spatial disorder into a random field acting on the Ising order, 
driving the transition into the paramagnetic QH phase. Fi- 
nally, the same coupling permits strain to act as a valley Zee- 



man field, and anisotropy to serve as a probe of transport — 
both important to experimental studies of nematic ordering. 

This perhaps a good place to observe that other situations 
in which the valley ordering involves higher symmetry have 
been studied in the past. A good example is graphene^^ : here, 
the Dirac dispersion is identical and to good approximation 
isotropic in the two valleys, and thus the em erge nt symme- 
try is SU{2). Another case of historical interesPElis the (1 10) 
surface of Si, the original example of valley QHFM. Although 
in bulk Si the valleys indeed have substantial anisotropy ori- 
ented along different axes, the two valleys that survive in the 
low-energy dispersion upon projection into the (110) plane 
have identical anisotropics; therefore, the symmetry here is 
again SU(2). In these cases, quite different phenomena 
emerge, such as low-energy skyrmionic 'valley textures' and 
gapless neutral Goldstone modes associated with the breaking 
of the continuous valley pseudospin symmetry. Furthermore, 
the equivalence of the anisotropy in the two valleys strongly 
diminishes the role of disorder, as it can no longer serve as a 
valley- selective random field; it thus does not couple directly 
to the pseudospin index and can act upon it only via the charge 
sector. 

Returning to our central topic, it is clearly desirable to find 
other instances of Ising-type valley QHFMs. In closing, we 
would like to flag a few examples as worthy of further study. 
The first of these. Si (111) heterostructures^^ ^^, possess six 
inequivalent valleys; these split into three pairs, with the dis- 
persion in the two valleys belonging to each pair exhibiting 
identical anisotropy. Considerations analogous to those pre- 
sented above suggest that the resulting QHFM should have 
an SU{2) x Zs symmetry, where the SU (2) rotates between 
the two valleys within a pair, and the discrete Z3 index acts 
between the three pairs, once again intertwining spatial and 
pseudospin rotations. While this is a considerably more intri- 
cate symmetry structure than the one considered in this paper, 
a naive expectation is that now Ising ordering occurs at z^ = 2, 
and involves filling the lowest Landau level in both valleys be- 
longing to a pair thus breaking the Z3 symmetry. While trans- 
port measurements suggest that u = 2 does indeed exhibit pe- 
culiar behavior, Z3 symmetry breaking is more subtle from the 
point of view of application of strain and the resulting trans- 
port anisotropy. Fully characterizing the symmetry-breaking 
transition, its experimental consequences, and potential exper- 
imental probes, remains an open question. A second exam- 
ple, bilayer graphene, would at first sight appear to exhibit the 
SU(2) symmetry of its monolayer cousin; however, the in- 
clusion of 'trigonal warping' effects into the band structure^S 
could break this down to an Ising symmetry. Once again, we 
defer a detailed study of this to future work. Finally, a far 
more speculative example is the possibility of similar transi- 
tions occuring in low-carrier-density systems in three dimen- 
sions; recently, transport experiments in BismutlP^M^ have 
demonstrated orientational symmetry breaking in the pres- 
ence of a magnetic field that is not too far from the quantum 
limit, which could be consistent with some valley-ordering 



scenarios.^ 
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